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Abstract. First, an abstract scheme of constructing biorthogonal rational 
systems related to some interpolation problems is proposed. We also present 
a modification of the famous step-by-step process of solving the Nevanlinna- 
Pick problems for Nevanlinna functions. The process in question gives rise 
to three-term recurrence relations with coefficients depending on the spectral 
parameter. These relations can be rewritten in the matrix form by means of 
two Jacobi matrices. As a result, a convergence theorem for multipoint Pade 
approximants to Nevanlinna functions is proved. 



1. Introduction 

Moment problems as well as interpolation problems arise in a wide range of 
mathematical and physical sciences (see [I], [5], [21] )■ The classical power moment 
problem can be formulated as follows. 

The Hamburger moment problem. Given is a sequence {cj}j°. of real num- 
bers. Find a positive Borel measure da on R such that 

t j da{t), j = 0,1,.... 

In a view of the Hamburger-Nevanlinna theorem (see [1]), the moment problem 
is equivalent to the problem of finding the Nevanlinna function ip(X) (= J R 
having the following asymptotic expansions 

< ^( A ) = -y- A l a Stt + °( a ^Tt) (A = iy,y^+oo) 

for all n G Z + := MU {0}. The latter problem can be solved by means of the Schur 
algorithm. This algorithm leads to the J-fraction 
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A - a 



b 2 



where aj are real numbers, bj are positive numbers [T] (see also [2]). Recall that the 
theory of J-fractions is also related to the Pade approximation theory and the theory 
of orthogonal polynomials. Under some natural conditions it is possible to say that 
all these theories (J-fractions, Pade approximation and orthogonal polynomials) are 
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equivalent to one another. On the other hand, the J- fraction (II. ip generates the 
following infinite Jacobi matrix 

/ a b Q \ 
6 ai h 

3 = h 

In fact, the Jacobi matrix J' is a key tool for analyzing the moment problem as 
well as the Nevanlinna function if via operator methods. For example, using Jacobi 
matrices techniques one can prove convergence results for Pade approximants to 
Nevanlinna functions (see, for instance, [5] , [57] ; see also [TT] where generalized Ja- 
cobi matrices associated to indefinite moment problems for generalized Nevanlinna 
functions are considered). 

The main goal of the present paper is to generalize the above-described scheme 
to the case of Nevanlinna-Pick problems in the class of Nevanlinna functions. Our 
approach is based on the relations between the theory of multipoint Pade approxi- 
mants (Pade interpolants), the theory of biorthogonal rational functions, and gen- 
eralized eigenvalue problem for two Jacobi matrices [31], [28], [29], [32] . 

In theory of biorthogonal rational functions the so-called continued fractions of 
Rn type appear. These continued fractions were introduced and studied by Ismail 
and Masson [T7] . Nevertheless, note that continued fractions of the same type were 
considered earlier in problems connected with rational interpolation problems (see, 
e.g. [30], [10]). It appears that the continued fractions of the Rn type are closely 
connected with the diagonal Pade interpolation problem from one side and with 
the theory of generalized eigenvalue problem for two Jacobi matrices on the other 
side [21] . In turn, both problem are equivalent (under some natural conditions) to 
theory of the biorthogonal rational functions (BRF) [33] , [32] • Note that theory of 
orthogonal rational functions studied and developed in [9] can be considered as a 
special case of theory of BRF (for details see, e.g. [32]). 

2. Pade interpolation and biorthogonality 

In this section we present basic facts concerning Pade interpolation and corre- 
sponding biorthogonal rational functions. We follow mostly [31], [32] but some of 
the result appear to be new. 

Let monic polynomials P n (z) = z n + 0(z n ^ 1 ) satisfy the Rn type recurrence 
relation 

(2.1) P n +i(z) + (a n z + f3 n )P n {z) + r n (z - a n ){z - b n )P n -i{z) = 
with initial conditions 

P = 1, Pi = z - A) 

Monic property of the polynomials P n {z) assumes the restriction upon the recur- 
rence coefficients 

(2.2) a = -1, a n + r n + 1 = 0, n = 1, 2, . . . 

In what follows we will assume that r n ^ 0, n = 1, 2, . . . (nondegeneracy). 
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Introduce the polynomials 

n n 

A = B = 1, A n (z) = Y[(z a k ), B n {z) = ]J(z b k ) 

k=l k=l 

As shown by Ismail and Masson [17j there exists a linear functional a defined 
on all rational functions (without a polynomial part) with the prescribed poles 
Oi, &i, a-ii &2j • • ■ by the moments 

(2.3) c nm = <?{ -* , } u rr l] n, to = 0,1,2, .. . 
such that the orthogonality relation 

(2.4) 4#?1t\}= > i = 0,l n-1, 

I A n (z)B n (z) ) 

holds, where qj(z) is any polynomial of degree not exceeding j and 

l A„(z)B„(z) J 

The normalization coefficients K n satisfy the recurrence relation [17j 

(2.5) Kn+i + a n K n + r n n n -i = 

It is important to note that, in contrast to the case of the ordinary orthogonal 
polynomials, we can take two first coefficients kq, K\ as arbitrary parameters. Then 
all further coefficients k 2 , k 3 , . . . are determined uniquely through (|2.5[) . 

Note also that if for some n — n > 1 we have K no — K no -i then from (|2.5p and 
(|2.2[) it follows that n na +i — K n a = K no -i and hence we then have K n = n no for all 
n > no — 1. Moreover, we also have from (|2.5[) 

r»o-l( K TlD-l _ K "o-2) = 

Due to our assumption r n ^ we have k„ _2 = n na -i — K n - Repeating this 
process we arrive at condition 

(2.6) «i = Kq. 
We thus have 

Proposition 2.1. Condition K no — K no —i for some no > 1 is equivalent to the 
condition kq — n±. In this case we have K n = const for all n = 0, 1, 2, . . . 

This case will be considered as a degeneration and in what follows we will assume 
that K\ 7^ kq. Then from this proposition it follows n n ^ « n _i for n — 2,3,.... 
Moreover we will assume that K n ^ for all n = 0, 1, 2, 

Introduce the rational functions [31j 

(2J) Hn {Z) A n ( Z y Kn [Z) B n {z) 

It is assumed that zeroes of polynomials P n (z) do not coincide with points ai,bj, 
so rational functions Rn\z) and Rn\z) have the [n/n] type. Rational functions 
R n (z)^ have prescribed poles ai, a2, . . . , a n and rational functions Rn\z) have 
prescribed poles b\, b 2 , . ■ . , b n . 

These functions satisfy obvious recurrence relations 

(2.8) (z - a n+l )R^ ] +1 {z) + (a n z + f3 n )R^\z) + r n (z - b^R^z) = 
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and 

(2.9) (z - b n+1 )R^ +1 {z) + (a n z + (3 n )R%\z) +r n (z- a^R^z) = 

On the other hand, these recurrence relations can be rewritten in terms of the 
generalized eigenvalue problem (GEVP) [31] 

,hR^ = zJ 2 R^ 

and 

J 3 R^ = zJ 2 R(V 

where R^ is an infinite-dimensional vector with components {Rq , R± , ■ ■ ■ } (as 
well as W- 2 ^ ) and J\ , J 2 , J3 are 3-diagonal ( Jacobi) matrices which entries are ob- 
vious from the above recurrence relations for Rn\ Rn^ ■ As was shown in [31] the 
GEVP leads naturally to theory of biorthogonal rational functions associated with 
the polynomials P n {z) of the Ru-type. Here we propose a more simple scheme of 
construction of the pair of biorthogonal rational functions. 

Introduce the rational functions U n (z) and V n (z) by the formulas: 

(2.10) U n {z) = R^{z)-^ n R ( t ] ) _ 1 {z), V n (z) = RW(z)-Z n R^ 1 (z) 

where = k„/k„_i (it assumed that £0 = so that Uq = Vq — 1). Clearly, the 
rational functions U n (z) have the poles a\, a 2 , ■ ■ ■ , a n and the rational functions 
V n (z) have the poles bi,b 2 , 
We have 

Theorem 2.2. The rational functions \2.10\) form a biorthogonal system with re- 
spect to the functional a: 

(2.11) o-{u n (z)V m (z)} = h n 5 nm , n,m = 0, 1, ... 
where the normalization coefficients are 

hn = ( K n— 1 ~ K n) 

The proof of this theorem is direct by using orthogonality relations (|2.4p . 

Note that the normalization coefficient is nonzero h n ^ due to our assumptions 
on nondegeneracy kq ^ n\ and k„ ^ 0. 

We can give an equivalent definition of the functions U n {z) and V n {z) using the 
determinant expressions: 

Cn,0 
Cn,l 



Cro,n— 1 

Co,n 
Cl, n 



(2.12) 



U n (z) 



Pn{a n ) 



coo 

Coi 



ClO 

en 



C0,n-1 Ci.„_i 

1 Aj\z) 



(2.13) 



V n (z) 



Coo 

ClO 



coi 

Cll 



Cn-1,0 Cn-i^ 

1 B^(z) 
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where 



(2.14) 



A„ = 



coo 

ClO 



COI 



Co,n-l 
Cl,n-1 



1,0 ^n— 1,1 • ■ ■ ^-n— l.n— 1 

(It is assumed that Ao = 1)- In what follows we will assume that 

A„^0, n= 1,2,3,... 

(this is another nondegeneracy condition). 

Formulas (|2.12|1 , (|2. 13|) follow directly from definition of moments (|2.3j) . In order 
to obtain appropriate coefficients in front of determinantal expressions (|2.12|) , (|2.13p 
we can present expression for the rational function U n (z) in the following form 



U n (z) 



\ -v 7nfc 



^ A k (z) 

k=0 y ' 



The leading term in this sum is 

Inn = U n (z)A n (z) \ z=an 

On the other hand we have from the explicit expression (|2. 10|) 

U n (z)A n (z) \z= an =P„(a„) 

whence we obtain the factor P "^°"^ in front of determinant of the formula (|2. 12[) . 
Similarly we obtain the factor P '^ b "- ) in front of determinant of the formula (|2.13[) . 

Note also that from the determinantal formulas (|2.12| . (]2 . 13[) it follows directly 
that 

(j{Un(z)V m (z)} = P n (a n )Pn(bn) $nm 

Comparing with (|2 . we obtain an interesting relation 

A, 



(2.15) 



Pn(a n )P n (b n J 



K n(l K n/ K n— l) 



A n +i 

From this relation it follows that condition 
(2.16) Pn{a n )Pn{b n ) + 

guarantees nondegeneracy conditions n n ^ 0, k„ ^ K n -i and A„ 7^ 0. Thus we 
will assume that condition (|2.16|) holds. It is instructive to consider what happens 
when condition (|2. 16|) doesn't hold. Assume e.g. that P n (a n ) = for some n. 
Then the rational function Rn\z) = P n (z) / A n (z) has the order [n — 1/n — 1], i.e. 
it has poles a\, 0,2, ... , a n -i. Corresponding rational function U n (z) constructed by 
(12.10)) will also have poles a\, a^, ■ ■ . , a„_i which means a degeneration (absence of 
the pole a n ). 

We can present rational functions U n (z) and V n (z) in the form 



(2.17) 



U n (z) = 



S n (z) 



(i-t n )A n ( z y 



V n (z) = 



T n (z) 



(i-e n )B n (z)' 
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where S n (z) = z n + 0(z n 1 ) and T n (z) = z n + 0(z n 1 ) are monic polynomials of 
degree n. Polynomials S n (z),T n (z) are expressed in terms of polynomials P n (z): 

(2.18) 

Of) - ~ £n( Z - a n )Pn-l{z) . , = Pn(z) - £,n(z - 6 7I )P ra _i(z) 

1 sn 1 C,n 

Moreover So=T = l. 
We have 

Proposition 2.3. Polynomials S n (z),T n (z) satisfy a system of first- order recur- 
rence relations 

S n +i(z) = vi 1] {z - b n )S n {z) + v£\z - a n )T n (z) 

(2.19) T n+1 {z) = v^(z - b n )S n {z) + v&(z - a n )T n (z), n = 1,2, . .. 

where 







£,nCn+ l&n+l Tn&n 




r n (b n - a n ) 


i/ 2 ) 


Cnftn 






r n (a n - b n ) 


u m 


£,nPn 


— £n£n+l6 n +l — r nO"n 




r n (b n - a„) 










r n {a n - b n ) 


Vrf* +Vn^ = 1 which is necessary for 



to be monic. 

Proof. It is sufficient to write down 

(9 0n\ TT f~\ - j3 »+l( Z ) c P n( z ) _ S n+1 (z) 
(Z.ZV) U n+ i[Z) — — — — Cn+l- 



A n+ i(z) A n (z) (1 - £ n +i)A n+ i 

and express P„+i(z) in terms of P n (z) and P n -i(z) using recurrence relation 
(12. ip . This allows one to obtain an expression of P n (z) in terms of polynomials 

S„(z),S„ + i(z): 

(2.21) P n (z) = C } (S n+ i(z) -(z- b n )S n {z)) , 



where 



£(1) _ r «(l - £n) 



Analogously 

(2.22) P n (z) = Ci 2) (T„ +1 (z) - (z - a n )T n (z)) , 



where 



c 



(2) _ r n(l ~ £«) 



r n (a n - 671+1) - + a„6 n+ i) 

Then we arrive at relations (|2.19p . 

Note that relations (|2~T9f (as well as (|2~2Tj) . (|2^22]) ) do not hold for n = because 
coefficients Uq as well as ao,&o are not defined. However, relations (|2. 19|) will be 
valid for n = if we put 

So = To = 1 
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and 

(2.23) Co = 

K - Kl 

whereas ao, fro & n d fo can be arbitrary parameters. Indeed it is elementary verified 
that in this case we have from relations (|2.19[) for n = 

(2.24) &i(z)=z-\ , Ti(z)=z-\ 

K ~ Ki K — K l 

which is compatible with expression for S'i(z),Ti(z) obtained from (|2.18|) for n = 
1. Thus we can assume that relations (j2.19|) are valid for n = 0, 1,2, . . . under 
condition (|2.23[) . Note that this condition is formally equivalent to condition 

k - Kl 

K-i = 

ro 

obtained from recurrence relation l]2.5[) if one puts n = (with arbitrary nonzero 
ro). Equivalently, we can assume that for n — coefficients Un take the values 
w (i) _ kq(/3q - gp) + Aci(a - oi) ( 2 ) _ k (/3 - 6 ) + «i(fr - gi) 



and 



° (&o - ao)(«0 - Kl) ' ° (OO - &o)(«0 

(3) _ ^o(A) - go) + Kijag - h) (4) _ k (/3q - 6 ) + Ki(b - bj) 



{b - a )(K - ki) !! (a - &o)(ko - Ki) 

Vice versa, one can show that starting from the system (|2.19p with b n ^ a n , n = 
0,1,... and with initial conditions To = So = 1 one construct a pair of biorthogonal 
functions U n (z), V n (z) by (j2~T7|) [29] . 

The Pade interpolation problem [6] (sometimes called the Cauchy-Jacobi, Newton- 
Pade or multi-point Pade approximation problem [24], [5]) consists in finding a pair 
of polynomials P n {z) 1 Q m (z) such that 

(2.25) Y s P n (z s ) = Q m (z 3 ), s = 0,1, 2,..., n + m, 

where Y s and z s are two given complex sequences (s = 0,1,2,...). The degrees 
of polynomials P n {z), Q m {z) are less or equal to n and m correspondingly. The 
rational function 

/ s _ Qm{z) 

rmn[Z) ~ P n (z) 

is called the rational interpolant corresponding to the sequences Y s and z s . 

We will consider only the so-called normal case of the Pade interpolation problem 
[H] meaning that the degrees of polynomials P n (z), Q m (z) are exactly n and m and 
there are no common zeros of polynomials P n (z) and Q m (z). In the normal case 
we have for every pair (n, m) the conditions [6] 

In practice, it is assumed that Y s = F(z s ) for some given function F(z) of the 
complex argument z. In this case formula (|2.25p gives a rational interpolant r mn {z) 
of the function F(z) for the given sequence z s of interpolation points. Note that 
when all interpolation points coincide z s — zo, s = 0, 1, 2, ... , then we have the 
ordinary Pade approximation in the point zq. The set r mn {z), m,n = 0, 1, 2, . . . is 
called the Pade interpolation table for the function F(z). 



8 



MAXIM S. DEREVYAGIN AND ALEXEI S. ZHEDANOV 



Consider the so-called diagonal string [29], [32] in the Pade interpolation table, 
i.e. the set r n -i in (z), n = 1,2, . . . . This means that we are seeking a solution of 
the problem 

(2.26) F (z s )= Qn -] {z ;\ S = 0,l,2,...2n-1 

Pade interpolants for the diagonal string satisfy simple orthogonality properties 

(2.27) [z ,z 1 ,...z 2n _ 1 ]{zJ> n (z)} =0, j = 0,l,...,n-l, 
where 



(2.28) ...,*„_,](/«} ^ rK _, 0)({ 

is the divided difference of the order 2n—l from the function f(z). It is assumed that 
the integration contour T avoids all singularity points of the function f(z). Note 
that formula (|2.28[1 is called the Hermite form of the divided difference operation 
0. 

Orthogonality relation (|2 . 2T[) can be extended to biorthogonality relation for two 
rational functions U n (z),V n (z) as follows. Consider the diagonal Pade interpola- 
tion problem for the same function F(z) but with slightly modified interpolation 
sequence 

(2.29) F(z s ) = 8 = o, i, 2 , . . . , 2n - 2, 2n 

Pn{z s ) 

(i.e. for the given n we have 2n interpolation points as in the previous scheme (|2.26[) . 
but the final point Z2 n -i is replaced by Z2 n )- Construct the rational functions 
(2.30) 

T r ( x „ , , Pn(z) 

U„(Z) = r r, V n (z) - 



(z - zi)(z - z 3 ) . . . (z - z 2n -i) ' (z - z 2 )(z - Zi) . . . (z - z 2n ) 

Then the biorthogonality relation 

(2.31) [z ,zi, . . . z 2n -i\ < >=h n d nm , n,m = 0,1,2, ... 

I z - z ) 

holds with some normalization constant h n ^ [29], [32]. It is easily verified that 
polynomials P n (z) and P n {z) satisfy the Rji recurrence relations (|2.1[) whereas 
the rational functions U n (z),V n (z) satisfy the generalized eigenvalue problem of 
type (|2.8[) . Thus the generalized eigenvalue problem for two Jacobi matrices is 
related with the diagonal Pade interpolation problem. For further development 
and generalizations of this subject see [29], [32], [23]. 

3. Nevanlinna-Pick problems 

In this section we propose a modification of the famous step-by-step process of 
solving the Nevanlinna-Pick problem in the class of Nevanlinna functions p] , [2] . 

First, let us recall that a Nevanlinna function is a function which is holomorphic 
in the open upper half plane C+ and has a nonnegative imaginary part in C+. Let 
N[a,/3] denote a class of all functions <p having the representation 

r p da(t) 



(3.1) ip(X) 
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where da(t) is a finite measure. A function of the class N[a, f3] is called a Markov 
function. Clearly, a Markov function is also a Nevanlinna function. Moreover, if the 
singularities of the Nevanlinna function ip are contained in [a, (3] then if G N[a,/?] 
(see, for instance, [1]). Let us consider the following Nevanlinna-Pick problem. 
Problem NP[q,/5]. Given are two infinite sequences {zk}% x ' = o, { w k}kLo ( z k £ C+). 
Find a function <p £ N[a, f3] such that 

ip{z k ) = w k , k = 0,1,2, 

As is known (see [2]), the problem NP[a,/3] is solvable if and only if the Her- 
mitian forms 

( 32 ) J2 w j( z J ~ a )- Wkjzk -«) ,^ J2 Wj{f3 - z 3 ) -w k (P -z k ) 
i,fc=o Zj ~ Zk ' j,k=o Zj ~ Zk 

are nonnegative definite for all N £ Z + . 

It is also natural to consider the truncated Nevanlinna-Pick problem. 
Problem NP([a,/3],n). Given are two finite sequences {2k}jJ_ , {wk}^ =0 (z k £ 
C + ). Describe all functions <p £ N[a,/3] satisfying the property 

ip(z k ) = w k , k = 0, 1, . . . ,n. 

Note that the problem NP([a,/?],n) is solvable if and only if the Hermitian 



forms (|3.2p are nonnegative definite for N = 0, 1, . . . , n. 

The algorithm of solving the Nevanlinna-Pick problems in question is based on 
the subsequent statement. 

Lemma 3.1. Let tp G N[a, j3] and let z £ C+ be a fixed number. Then there exist 
numbers E K and b > such that the function r defined by the equality 

(3 - 3) ^ (A) = ~ a (2) A - a(D + 6 2 (A - z)(X - z)r(A) 

&eZon 5 s to N°[a,,3] := N[a,/3] U {0}. 

Proof. Setting $(A) := --— , define the function 

(3.4) *(A) = : *Z±. 

$(A) -$(z) A-z 

Due to the Schwartz lemma, we have that 

|*(A)| < 1, ImA>0. 
So, the function defined from the relation 

(3.5) = 

#l(A) +2 

is a Nevanlinna function. Plugging (|3.5[) into (|3.4[) . one obtains 



(3 6) $ (A) = ^( A )( 2A ~ g - g) ~ - z ) - gCgX* - f) 

$(A)(z-z)-$(z)(A-z) + $(z)(A-z) 

Now, let us consider the following function 

T / \ \ t ,,s . 2A - (z + z) ~ 2A-(z + z) 

i(A) := *i(A) + i = *i(A) + — i '-. 

z — z 2 1m z 
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Obviously, is a Nevanlinna function. Taking into account (13. 6j) . admits the 
following representation 

Im$(z) (A-z)(A-z) 



(3.7) *!(A) 



Imz $fV) Im$ ( z ) \ I IslWj 



Im js Im z 

Finally, introducing 

1 , T , Im$(z) r9 N Im$(z) m m Im<&(z)z 

t(A) = — - — -— — G N, 6 = V ; > 0, = v ' e R, a« = — L_L e 

W^A) Imz Imz Imz 



one can easily transform (|3.7|) into (|3 . 3|) . To complete the proof, it is sufficient to 
observe that, due to (|3.1[) and (|3.3[) . all singularities of r are contained in [a, /3]. □ 



Remark 3.2. The transformation (|3.3[) could be viewed as a substitute for the 
Schwartz lemma. A similar to (|3.3[) transformation for Caratheodory functions was 
proposed in [ID] . 

Remark 3.3. Substituting A for z and z in (|3.3p we get 

= — f2T~ — riY> ^ = — TTP" — fir 

Expressing from the above relations and a^ 2 \ one can obtain the following 
formulas 

, , Im —t^t , , Im —t^ i 
(3.8) a& = a« = ^z + 



Imz Imz V 3 ! 2 ) 

It is easy to see that the numbers a*- 2 -* are uniquely determined by 

Further, equality (|3.3|) can be rewritten as follows 



i +0 (2) A _ (i) 

^ ^( a > = - *Ta-«)(a-*) ■ 

In fact, the number 6 can be chosen arbitrary. So, to be definite we always choose 
b > in the following way 

P 



b 2 = da(t) 

J a 



In this case, the function t possesses the integral representation (13.11) with a prob- 
ability measure. 

Remark 3.4. It also easily follows from the theory of generalized Nevanlinna 
functions (see [12], [13], [20]) that the right-hand side of (|3.9p is a Nevanlinna 
function. 

Remark 3.5. By comparing the first terms in asymptotic expansions of the right 
hand side and left hand side of (13.31). we see that 




(3.10) a^> = / da(t) + b 2 . 



Now, we are in a position to solve the problem NP([a, 0\, n). Let the given prob- 
lem NP([a, /3],n) be solvable and let ip be a solution of the problem NP([a, 0\, n). 
Due to Lemma 13.11 tpo := <p admits the following representation 

(3.11) ip(X) = -—^ 



ag'X-ctf' + b 2 (X-z )(X-z a ) Vl (X) 
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where pi € N°[a, /3]. From (I3.11|) we see that 

!Pi(A) 



1 

RAJ 



a A 



6g(A - z )(A - zq) 



So, if ipi ^ then it is a solution of the problem NP([a, 0], n— 1) with the sequences 
{zfc}£=i and {^i^lLi' where 



w 



(i) 



^i(zjfe) 



— + aiP Zk 



^(^fc _ ^o)(2fe - «o) 



Therefore, the original problem NP([a, (3], n) is reduced to the problem NP([a, 0\, n 

1) . Similarly, the problem NP([a, (3], n— 1) can be reduced to the problem NP([a, 0), n— 

2) and so on. Finally, one has a sequence of the linear fractional transformations 

^( A ) = — w- 



■ W +b^(X-z J )(X-z J )^ 
having the following matrix representations 



a- >X-a) 



l(A) 



(i = o, i, . 



(3.12) 



Wi(A) 







6j(A-Zj) 



1 



bAX-Zj) 
o (2 "a 



(j = 0,l,...,7l). 



bj(X-Zj) J 

If the above-described algorithm consists of exactly n + 1 steps then we say that the 
problem NP([a, /?],n) is nondegenerate. So, we have proved the following theorem 
which gives the complete solution of the problem NP([a, 0\, n). 

Theorem 3.6 ([lj). Any solution ip of the nondegenerate problem NP([a,/3],n) 
admits the following representation 

w n (X)r(X) + w 12 (X) 



(3.13) 

where r G N [a, (3\ and 
(3.14) W [0 ,„](A) 



p(X) = 



W2l(X)r(X) + W 2 2(X) 



K'(A))?,- =1 := W (A)Wi(A) . ..W„(A). 



It should be also remarked that Wj is the Blaschke-Potapov factor [3] , [26] . 

4. Rational systems related to Nevanlinna-Pick problems 

Let us suppose that the given Markov function has the integral representa- 
tion (|3.1[) with a probability measure da which support contains infinitely many 
points, i.e. 

da(t) = 1. 



Let us also suppose that for the given sequence {zk\^ =0 there exists 5 > such 
that 



(4.1) 



lmz k >S, k = 0, 1,2, . 



In this case, it follows from the uniqueness theorem for analytic functions that the 
numbers Zk and Wk '■= p>{zk) (k € Z+) uniquely determine the function ip. So, the 
Nevanlinna-Pick problem with the data {zk)kLoi i w k}kLo nas a un iQ ue solution. 
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Since tp is not rational the given data give rise to infinitely many steps of the 
step-by-step process. So, we have infinitely many linear fractional transformations 
of the form (|3.3|) which lead to the following continued fraction 

(4.2) - ' 



(2) , _ (1) fr5( A - z o)( X - z o) 



a y ,'X-aY 

6g(A- z )(X-z ) bfiX-zx^X-zx) 



„(2) A - J2K (1) - J2K (1) • 

Aq a — &q A — a^ Q>2 A — 

The continued fraction (|4.2p is an i?//-fraction (see [T7]). Consider the (n + l)-th 
convergent of the continued fraction (|4.2|) 

1 &g(A - z )(X - z ) ^-i( A - 2„-i)(A - z„_i) 



Aq A — CLq A — &n A — (In 



It is obvious that R n is a solution of the problem NP([a, f3], n), i.e. the following 
equality holds true 

Rn(zk) = w k = <p(zk), k = 0, l,...,n. 

Definition 4.1. The [L/M] multipoint Pade approximant for a function tp at the 
points {ak}kLi is defined as a ratio 

1 { ' BIVMI(A) 

of two polynomials AS L ' M \ B^ L I M ^ of formal degree L and M, respectively, such 
that 

f [L/M] {a k ) = tp(a k ), /-• 1 /. • M I I. 

Since R n is real, the rational function R n is the [n/n] multipoint Pade approxi- 
mant for tp at the points {oo, zq, zq, . . . , z n ,"z„, ■ • ■ }■ 

It is well-known that to every continued fraction there corresponds a recurrence 
relation. In particular, for the continued fraction (|4.2[) a recurrence relation takes 
the following form 

(4.3) u j+1 - [afx - af^uj + b^X - z J _ 1 )(A - z^uj-i =0 (j£ N). 



Define polynomials of the first kind Pj (A) as solutions Uj — Pj (A) of the system (|4.3[) 
with the initial conditions 

(4.4) 1*0 = 1, u 1 = ctf ) \-c$ ) . 

Similarly, the polynomials of the second kind Qj(X) are defined as solutions Uj = 
Qj{X) of the system (I4.3|) subject to the following initial conditions 

(4.5) u Q = 0, «i = -1. 

Note that in our setting (|2.4I) is transformed into the following orthogonality 
relations (see also [15] ) 
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It follows from the theory of continued fractions that Rn(X) = ^ ( se e, for 

details, [18]). 

The recurrence relation (|4.3[) can be renormalized to the following one 



(4.7) - X)u j+1 - [afx - af } )u 3 + bj-ify-i - A)%_! =0 (j G N), 



where 



6 . . . 6j-i(zo - A) . . . (zj_i - A) 



0' e N). 



Relation (|4. 6|) implies that 
(4.8) 
Now, setting 

' rP V +1 da{t) 



£o = 0, & 



" - 1 
P*+i (t) — —d(i{t) = 0, j = 0, . . . , n. 

£ - 2, 



t j da{t) 



a \t-z \*...\t- Zj \* 



a \t-Z \ ...\t- Zj-l\ 



(j G N) 



one can see that the simple linear combinations Pj — £jP,-_i (j G Z + ) give orthog- 
onalization of the system 



A — zo A — z\ 

of rational functions (see Theorem 12.21 see also [7])- It should be also noted here 
that systems of orthogonal rational functions related to Nevanlinna-Pick problems 
were proposed in [8], [22], [25] (see also [9]). 

Further, relation (|4.7p can be rewritten as follows 



(4.9) ZjbjUj+i +ay'v,j + Zj-ibj-iUj-i = X(bjU J+ i +a[ 2) Uj +b j _iu J ^ 1 ) (j G N) 



The system (|4.9[) gives us the possibility to rewrite the Cauchy problem (|4.3|) . (|4.4|) 
in the matrix form 



4 oo)^ A ) = Mo/oo^)- 



where ^(A) = ( P (A), Pi(A), . . . , PAX), . 



J, 



(i) _ 



/ ajp z b 

Zo b a{' zih 
h (!) 

V 



J, 



and 



(2) _ 



• "0 



[0,00) 



b af ] bi 

h ( 2 ) 

61 a 2 



We denote by ^ n i the space of (n + 1) vectors with the usual inner product. Define 
a standard basis in £^ n , by the equalities 

e j = {^,fc}Lo> j = 0, 1,...,ti. 
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^\jM ^UM> ^ e a submatrix of J\n t00 \ (JrgL)) corresponding to the linear 
subspace spanned by the vectors e;, . . . , (0 < j < k < n), that is, 



V 



A4 2) h o\ 



j, 



r(2) 



(2) 

b'.fc] 



V 



,(2. 



Proposition 4.2. TTie matrix Jj n j is positive definite for all n G Z + . 
Proof. Let us consider the Hermitian form 

(4.10) (jfoij^e) = 4 2) N 2 + 6ofo6 + &o&?i + a^\^\ 2 + ■■■ + a^ n \ 2 . 

Due to (|3.10p and our assumptions, we have that a!"p = 1 Therefore, one can 
rewrite the form (14. 10[) in the following manner 

= ICO I' + l&0£0 + 6 I 2 + • • ■ + |&n-l£n-l + in? + KU" > 0- 

(2) 

Thus, JL n j is a positive definite matrix. □ 

Finally, we should note that for the matrix JL ' n the following factorization 
holds true 

I af b \ ( 1 6o 

1 h 



(2) 

[0,oo) 



bo ai 2 bi 

h (2) 

bi a\ 



\ ( 1 

6o 1 



1 



1 V 



6i 1 



/ V 



5. m-FUNCTIONS OF LINEAR PENCILS 

In this section we give a matrix representation of multipoint Pade approximants 
for Markov functions. 

Definition 5.1. The function 

(5-1) ^(^-((^-Ajgi,)- 1 ^,^) 
will be called the m-function of the linear pencil Jy n ] ~ ^^[jli]' 

To see the correctness of the above definition it is sufficient to rewrite (|5.ip in 
the following form 



(5.2) 



From (15.21) one can conclude that 



is a Nevanlinna function. 



Proposition 5.2. T/ie m-functions TOu „i and m^ + i n ] are related by the equality 

1 



(5.3) 



l [iM 
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Proof. Consider the following block representation of the matrix Ju n ] — ^uli] 

t(1) _ X t(2) _ 

where B = (bj(~Zj — A),0, ...,0). According to the Frobenius formula pTHl Sec- 
tion 0.7.3] the matrix (J^L — AJ^) -1 has the following block representation 
(5.4) 

Plugging (|5T4|) into fBTT]) . one obtains (pT3j) . □ 
Corollary 5.3. T/ie following equalities hold true 

(5.5) m [0 ,»](A) = iZ„(A) = £ N[a,/?]. 

Proof. Relation (|5.3| implies that the rational functions m[o jn ] and R n have the 
same expansions into Ru -fractions. So, m[o, n ] = Rn- By using standard argu- 
mentation, from (|4.6p one can conclude that a// i/ie zeros of P n +i are contained in 
[a,/3] (see pQ, [E]). The latter means that the Nevanlinna function mm^i belongs 
toN[a,/?]. ' □ 

So, now one can say that R n is a solution of NP([a, /3], n). By using standard 
argumentation, from (|5.5|) we can conclude the following result. 

Corollary 5.4. The zeros of P n+ \ and Q n +i are interlace. 

Below, we will need the following statement. 

Corollary 5.5. The spectrum a ^^n] ("^pn]) 1 ) °f ^ e ma ^ r * :r ^[oijt^oi])" 1 * s 
contained in [a, (3] . 



Proof. From the formula for calculation of inverse matrices, (|5.1|) . and ()5.5p one 
can see that 

m ^(Jfij-Ajgij) Qn+1 (A) 



So, the statement immediately follows from Corollary 15.41 and the fact that all the 
zeros of P n +i are contained in [a, 0\. □ 

Remark 5.6. It should be remarked that, for the case of the Laurent orthogonal 
polynomials, a similar scheme with two matrices and m-functions were considered 

in a. 

6. A CONVERGENCE RESULT FOR MULTIPOINT PADE APPROXIMANTS 



The goal of this section is to prove an analog of Markov's convergence theorem 
by making use of the operator representation of multipoint Pade approximants. 
We begin with an auxiliary statement. 
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Lemma 6.1. The following inequalities hold true 

(( J [oi])" le o^o) < 1 (ne 
Proof. The proof is by induction. First, note that 



((^[n,n] 



_1 e„, e„) = = < 1 (n£ Z+). 



Suppose that ^(J^^ n j ) 1 ek+i 1 ek+ij < 1. It follows from the Riccati equation [TU 
formula (2.15)] (see also (|5.3p ) that 

(( J [fei]) lefc ' efe ) =_ (^ - / (2) — V 

a k - b l[( J [k+l,n])~ ^fe+l.efe+lj 



1 + ^ 



1 

< 1. 



□ 



Now, we are ready to prove the main result of this section. 



Theorem 6.2 ([15]). Let ip G N[a,/3] and let the sequence {zt}^° = i satisfy the 
condition (|4.1[1 . Then the sequence /["/"I = R n converges to tp locally uniformly in 
C\[a,/3]. 

Proof. We first recall the well-known estimate for the resolvent of self-adjoint op- 
erator J (for instance see [19j Theorem V.3.2]) 

Next, observe that the operator ^mLC^r^i) -1 is self-adjoint with respect to the 
following inner product 

(2) W„ -A „ ..cl (Tn+1 



(( j [m) lx >y) 



Taking into account the representation ()5.5|) . ()5.2p . the Cauchy-Schwartz inequal- 
ity, (|6.1|) . Corollary [531 an d Lemma [BTTJ we obtain 

(6 - 2) < ((^eo^ < 1 

~ dist(A, [a,0\) ~ dist(A, [a, [3])' 

It follows from (|6.2p and Montel's theorem that the family {R n } is precompact in 
the topology of locally uniform convergence in C \ [a, f3]. Note that 

R n {zk) = <p(zk), n>k. 
Thus, applying the Vitali theorem completes the proof. □ 



Remark 6.3. Theorem 16.21 was proved in [15] by means of another method. The 
rates of convergence of multipoint Pade approximants was also given in [15] . The 
operator interpretation of the rates of convergence and a more detailed analysis of 
the underlying linear pencil will be given in the forthcoming paper. 
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